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ABSTRACT 

We present results from multi-epoch spectral analysis of XMM-Newton and Chandra observations of 
the mini broad absorption line (BAL) quasar PG 1115+080. This is one of the few X-ray detected 
mini-BAL quasars to date that is bright enough in the X-ray band, mostly due to large gravitational- 
lensing magnifications, to allow in-depth spectral analysis. The present XMM-Newton observations of 
PG 1115+080 have provided the highest signal-to-noise X-ray spectra of a mini-BAL quasar obtained to 
date. By modeling the spectra of PG 1115+080 we have obtained constraints on the column density and 
ionization state of its outflowing absorbing gas. A comparison between these constraints over several 
epochs indicates significant variability in the properties of the outflowing absorbers in PG 1115+080. 
The depths of the high-energy broad absorption features in PG 1115+080 show a significant decrease 
between the first two observation epochs separated by a rest-frame timescale of ~ 1 year. This variability 
supports the intrinsic nature of these absorbers. Assuming the interpretation that the high-energy 
absorption features arise from highly ionized Fe XXV we constrain the fraction of the total bolometric 
energy released by quasars PG 1115+080 and APM 08279+5255 into the IGM in the form of kinetic 
energy to be e k = 0.64±° g (68% confidence), and ek = 0.09io!o5i respectively. According to recent 
theoretical studies this range of efficiencies is large enough to influence significantly the formation of the 
host galaxy and to regulate the growth of the central black hole. 

Subject headings: galaxies: active — quasars: absorption lines — quasars: individual (PG 1115+080) 
— quasars: individual (APM 08279+5225) — X-rays: galaxies — gravitational lensing 



1. INTRODUCTION 

In recent years there has been mounting evidence from 
both theoretical and observational studies for the impor- 
tance of quasar outflows in regulating the growth of su- 
permassive black holes, controlling the formation of host 
galaxies, and enriching the Intergalactic Medium (IGM). 
Models of structure formation in a A-CDM cosmology are 
not consistent with observations unless feedback, cither 
from star formation or active galactic nuclei, is included 
in the simulations. Recently, the potential importance of 
quasar outflows has been explicitly demonstrated in theo- 
retical models of structure formation and galaxy mergers 
that incorporate the effects of quasar outflows [e.g., Scan- 
napieco & Oh 2004 (SO04); Granato et al. 2004 (G04); 
Springel, Di Matteo, & Hernquist 2005 (SDH05); Hop- 
kins et al. 2005, 2006]. A basic assumption in the mod- 
els of SO04 is that all quasars host outflows. SO04 find 
that by choosing the fraction of the total bolometric en- 
ergy released over a quasar's lifetime into the Interstellar 
Medium (ISM) and IGM in the form of kinetic energy 
to be ek=0.05, they can successfully model the observed 
evolution of the 5-band quasar luminosity function be- 
tween redshifts of 0.25 and 6.25. In a recent study SDH05 
simulated the growth of black holes in gas-rich galaxies 
with and without the presence of accretion feedback. The 
feedback in their model is thought to occur through en- 
ergetic quasar outflows that interact with the gas of the 
host galaxy. It is assumed that about 5% of the radiated 



luminosity is thermally coupled via these outflows to the 
surrounding gas. They find that the growth of the black- 
hole mass is self-regulated and eventually saturates at a 
final value that depends on the initial amount of gas avail- 
able for accretion. These authors also found that feedback 
from stars and quasars in their simulations of galaxy merg- 
ers can heat and expel gas from the centers of the merged 
galaxies. This loss of gas during the final stages of the 
merger can halt nuclear starburst activity and cause the 
merger galaxy to become a gas-poor elliptical. 

With the advent of XMM-Newton and Chandra it has 
become possible to infer the kinematic and ionization 
properties of highly ionized X-ray absorbers in Seyfert 1 
galaxies (e.g., Kaspi et al. 2002; Kaastra et al. 2002; Net- 
zer et al. 2003). Approximately 60% of Seyfert 1 galaxies 
show outflowing X-ray and UV absorption by ionized gas 
with velocities up to f» 2,500 km s _1 (e.g., Crenshaw et 
al. 1999; Kriss 2002). The absorbing outflow properties 
of more luminous quasars appear to differ in several ways 
from those of Seyfert Is. The fraction of occurrence of 
X-ray and UV absorption has been reported to be lower 
in quasars (e.g., Ganguly et al. 2001; George et al. 2002); 
however, the claims for a lower occurrence of X-ray absorp- 
tion have been challenged by Porquet et al. (2004) who 
find that about half of their sample of 21 low-redshift PG 
quasars harbor ionized absorbers. UV spectroscopic ob- 
servations indicate that about 20% of quasars show broad 
absorption lines blueward of their resonant UV emission 
lines (e.g., Hewett & Foltz 2003). Outflow velocities of the 
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UV absorbers in BAL quasars have been found to be as 
high as w 60,000 km s" 1 . 

Recent X-ray observations of the mini-BAL quasar 
PG 1115+080 and BAL quasar APM 08279+5255 have 
suggested the presence of relativistic outflows of highly 
ionized (ionization parameter of the order of log£ = 3.5) 
absorbing material detected in the iron region (> 6 keV) 
with velocities of up to -0.4c (Chartas et al. 2002, 2003). 
The inferred hydrogen column densities ranging between 
2.g22-23 cm~ 2 and relativistic velocities of these outflowing 
X-ray absorbers imply mass-outflow rates that are compa- 
rable to the estimated accretion rates. The presence of 
massive, highly ionized, and high-velocity outflows from 
quasars indicates that these winds may be providing signif- 
icant feedback to the surrounding gas. Additional obser- 
vational evidence to support the presence of quasar feed- 
back came with the detection of high-velocity blueshiftcd 
absorption-line features in the X-ray spectra of several 
quasars and Narrow-Line Seyfert 1 galaxies (Reeves et al. 
2003; Pounds et al. 2003a, 2003b). We note, however, 
that some of these claims have been disputed in recent re- 
analyses of the data (e.g., Kaspi & Behar 2006; McKernan 
et al. 2004, 2005). 

The relatively small velocity spread observed for the UV 
absorption in PG 1115+080, compared to the typical range 
of 5,000 - 25,000 km s" 1 observed in BAL quasars, sug- 
gests that PG 1115+080 be classified as a mini-BALQSO 
(Turnshek 1988; Barlow, Hamann, & Sargent 1997). We 
will consider PG 1115+080 as a mini-BAL throughout this 
paper with the following model in mind. According to the 
"unification" model for BAL quasars most quasars have 
outflowing winds, however, BAL and mini-BAL quasars 
correspond to quasars with relatively large inclination an- 
gles. In particular, in BAL quasars it is commonly thought 
that our line of sight intersects a large portion of the out- 
flowing absorber whereas in mini-BAL quasars our line of 
sight intersects a shorter portion. According to this model 
mini-BAL quasars would be expected to have smaller ab- 
sorbing columns, and smaller velocity gradients than tra- 
ditional BAL quasars. 

In this work we present recent results from monitoring 
X-ray observations of the mini-BAL quasar PG 1115+080 
(z = 1.72). The goal of these observations was to monitor 
the time variability of the absorption features and thereby 
constrain the kinematic, ionization, and absorption prop- 
erties of the quasar outflows in this X-ray-bright mini-BAL 
quasar. Such constraints will allow us to estimate the con- 
tribution of quasar winds to the enrichment of the ISM 
and IGM and obtain a better understanding of the con- 
nection between black-hole and bulge growth in the host 
galaxy. 

Throughout this paper we adopt a A-dominated cosmol- 
ogy with H = 70 km s _1 Mpc" 1 , ft A = 0.7, and fi M = 
0.3. 

2. OBSERVATIONS AND DATA ANALYSIS 

PG 1115+080 was observed with XMM-Newton (Jansen 
et al. 2001) on 2001 Nov 25, 2004 June 10, and 2004 
June 26, for 62.9 ks, 81.2 ks and 86.3 ks, respectively. 
It was also observed with the Advanced CCD imaging 

4 The WWW sites listing the updates are located 
http : //xmm. vilspa. esa. es/external/xmm_sw_cal/calib/rel_notes/ 



Spectrometer (ACIS; Garmire et al. 2003) on board the 
Chandra X-ray Observatory (hereafter Chandra) on 2000 
June 2 and 2000 Nov 3 for 26.8 ks and 10.0 ks, respec- 
tively. The spectral analysis of the Chandra observations 
of PG 1115+080 has been presented in Gallagher et al. 
(2002) and Chartas et al. (2003). Because of recent signif- 
icant improvements in the calibration of the instruments 
on board Chandra and XMM-Newton since the publica- 
tion of thePG 1115+080 results, we have re-analyzed all 
observations. Updates on the calibration of Chandra and 
XMM-Newton are reported on the Chandra X-ray Cen- 
ter (CXC) and XMM-Newton Science Operations Centre 
(SOC) World Wide Web (WWW) sites, respectively]! We 
analyzed the XMM-Newton data for PG 1115+080 with 
the standard analysis software SAS version 6.1 provided 
by the XMM-Newton SOC. The Chandra observations of 
PG 1115+080 were analyzed using the standard software 
CI AO 3.2 provided by the CXC. A log of the observations 
that includes observation dates, observed count rates, to- 
tal exposure times, and observational identification num- 
bers is presented in Table 1. We note that the count 
rate and total number of counts (~31,430 counts from 
the three XMM-Newton observations) for PG 1115+080 
are the highest of any mini-BAL quasar X-ray spectrum 
observed to date. 

For the reduction of the XMM-Newton observations 
we filtered the PN (Striider et al. 2001) and MOS 
(Turner et al. 2001) data by selecting events correspond- 
ing to instrument PATTERNS in the 0-4 (single and dou- 
ble pixel events) and 0-12 ranges, respectively. Sev- 
eral moderate-amplitude background flares were present 
during the XMM-Newton observations of UM 425 and 
PG 1115+080. The PN and MOS data were filtered to 
exclude times when the full-field count rates exceeded 
20 cnts s _1 and 4 cnts s _1 , respectively. The extracted 
spectra from the PN and MOS were grouped to obtain a 
minimum of 100 and 40 counts, respectively, in each energy 
bin, allowing use of x 2 statistics. Background spectra for 
the PN and MOS detectors were extracted from source- 
free regions near UM 425 and PG 1115+080. The PN 
and MOS spectra were then fitted simultaneously with a 
variety of models employing XSPEC version 11.3 (Arnaud 
1996). The energy ranges used for fitting the PN and MOS 
data were 0.2-10 keV and 0.4-9 keV, respectively. 

For the reduction of the Chandra observations we used 
standard CXC threads to screen the data for status, grade, 
and time intervals of acceptable aspect solution and back- 
ground levels. The pointings of the observatory placed 
PG 1115+080 on the back-illuminated S3 chip of ACIS. 
To improve the spatial resolution we removed a ± 0"25 
randomization applied to the event positions in the CXC 
processing and employed a sub-pixel resolution technique 
developed by Tsunemi et al. (2001). 

In both the XMM-Newton and Chandra analyses we 
tested the sensitivity of our results to the selected back- 
ground and source-extraction regions by varying the loca- 
tions of the background regions and varying the sizes of 
the source-extraction regions. We did not find any signifi- 
cant change in the background-subtracted spectra. For all 
models of PG 1115+080 we included Galactic absorption 

at http://asc.harvard.edu/ciao/releasenotes/history.html and 
index . shtml , respectively. 
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due to neutral gas with a column density of Nn=3.5 x 

10 20 cm -2 (Stark et al. 1992). All quoted errors are at the 
90% confidence level unless mentioned otherwise with all 
parameters conservatively taken to be of interest except 
absolute normalization. 

2.1. Spectral Analysis of the XMM-Newton Observations 
of PG 1115+080. 

We first simultaneously fitted the PN, MOS1 and MOS2 
spectra of PG 1115+080 at each of the three epochs with 
a model consisting of a power law with neutral intrinsic 
absorption at z = 1.72 (model 1 of Table 2). These fits 
support the presence of an intrinsic absorber with col- 
umn densities of A H = 2.2±g"J x 10 21 cm" 2 , 4.6±g"J x 

10 21 cm" 2 , and 2.9±g;£ x 10 21 cm" 2 , for the three XMM- 
Newton observations of PG 1115+080. The fits are not 
acceptable in a statistical sense (see the reduced x 2 val- 
ues in Table 2). The fit residuals show significant absorp- 
tion at observed-frame energies of 2.5-5 keV. To illustrate 
the presence of these features and absorption features be- 
low 0.6 keV, we fit the spectra from observed-frame 2- 
2.5 keV and 5-10 keV with a power-law model (modified 
by Galactic absorption) and extrapolated this model to 
the energy ranges not fit (see Figure 1). For clarity we 
only show the higher S/N ratio PN data in Figure 1; how- 
ever, all fits were performed simultaneously using the PN 
and MOS1+2 data unless mentioned otherwise. The lower 
panels in Figure 1 show the Ax residuals between the best- 
fit power-law model and the PN data. We note that the 
apparent variability in the amplitude of the absorption 
features is not due to variations in the best-fit value of T 
for the spectral fits performed within the observed-frame 
ranges of 2-2.5 keV and 5-10 keV. The best-fit values of 
the photon indices for epochs 1, 2 and 3 were 1.91 jl ' 10 , 

1.88i^io; an d 1 -84lo.il- For the purpose of comparing the 
absorption residuals between epochs the photon indices for 
epochs 1,2, and 3 were set to T = 1.9. We proceed in fit- 
ting a variety of models to the data guided by the shape 
and location of these identified absorption residuals. As 
our first refinement we considered an absorbed power-law 
model with Gaussian absorption lines near the absorption 
features appearing between observed-frame energies of 2 
and 5 keV (see model 2 of Table 2). Hereafter, we refer to 
these absorption features as absl and abs2. The fits to the 
spectra of PG 1115+080 for epochs 1 and 2 contained two 
absorption lines, and the fits for epoch 3 contained one 
absorption line. The inclusion of these absorption lines 
in model 2 resulted in significant improvements of the fits 
compared to the previous ones in model 1 at the > 99.9%, 
99.2% and 93.6% confidence levels (according to the F- 
test) for epochs 1, 2 and 3, respectively. 

The remaining most significant contributions to the 
large values of \ 2 f° r model 2 arise from the residuals 
below 0.6 keV. To model these residuals we replaced the 
neutral absorber in our spectral model with an ionized 
intrinsic absorber (see model 3 of Table 2). These low- 
energy residuals are commonly detected in moderate S/N 
spectra of BAL quasars and are thought to arise from ab- 
sorption by multiple ions of O, Ne, Na, Mg, and Fe and/or 
partial covering. In particular, we used the absori model 
contained in XSPEC (Done et al. 1992). We note that the 
absori model is just a first approximation to what is likely 



a more complex situation. The temperature of the ionized 
absorber is not calculated self-consistently by proper ther- 
mal balance in the absori model but considered as an 
input parameter. We found a significant improvement in 
fit quality at the > 99.9% confidence level according to the 
F-test (see model 3 of Table 2). The best-fit rest-frame 
energies and widths of the absorption features for model 3 
for epochs 1, 2 and 3 are listed in Table 2. We note that 
the best-fit energies of these features are quite stable when 
trying different models for the low-energy spectral com- 
plexity. As an independent check of the accuracy of the 
fits that used the absori model we also repeated several 
of these fits using the warm absorber model XSTAR. XSTAR 
calculates the physical conditions and emission spectra of 
photoionized gases. In the current analysis we use a re- 
cent implementation of the XSTAR model that can be used 
within XSPEC. 

We found that the best-fit parameters of the low-energy 
ionized absorber derived using the XSTAR model were con- 
sistent with the values found using the absori model. We 
performed fits to the spectra of PG 1115+080 at the three 
epochs using both absori and XSTAR to compare results 
and performed the estimates of the confidence contours 
of the best-fit parameters using the more time-efficient 
absori model. 

We investigated whether the observed flux variability of 
the high-energy absorption features was consistent with a 
decrease in the strength of two absorption lines used to 
model the absorption features. The energies and widths 
of the absorption lines were held fixed between epochs. 
Specifically, we performed a simultaneous fit to all epochs 
incorporating model 2 of Table 2 with the difference that 
the energies and widths of the absorption lines were con- 
strained to be the same and only the normalizations of the 
absorption lines were allowed to vary. These fits were per- 
formed in the 1-10 keV band of the higher S/N PN data. 
We obtained a x 2 of 137.1 for 137 dof, and the best-fit rest- 
frame energies of the lines were 7.5 ± 0.3 keV and 10.0 ± 
0.4 keV. Allowing the energies and widths of the lines to 
vary resulted in a similar x 2 of 132.0 for 129 dof. We 
note that for the fits where the energies and widths were 
constrained to be the same between epochs we found that 
the best-fit values of the optical depths for the absorption 
lines at 7.5 keV were consistent with zero for epochs 2 and 
3. This was expected since an inspection of Figure 1 indi- 
cates that there are no residuals for epochs 2 and 3 near 
the rest-frame energy of 7.5 keV (left vertical dashed line) 
. We conclude that the strengths of the absorption lines 
for epochs 2 and 3 are too weak to infer any variability of 
their energies or widths. 

We next replaced the ionized absorber with a partial- 
covering model and repeated the fits for each epoch. We 
find that fits that include partial-covering (see model 4 of 
Table 2) provide x 2 values and quality-of-fit parameters 
comparable to those that included an ionized absorber. 
Models that included absorption edges for the high-energy 
absorption features of PG 1115+080 for epoch 1 were con- 
sidered and rejected in the current analysis, that includes 
the new calibration, as previously found in Chartas et al. 
(2003). The basic problem with such edge models is that, 
for plausible iron abundances, they predict too much ab- 
sorption at low energies. The inclusion of a broad (width 



4 



CHARTAS ET AL. 



allowed to vary in fits) Fe Ka emission line in all models 
in Table 2 for epoch 1 resulted in larger reduced \ 2 values 
compared to those with no line included; therefore, there 
is no evidence for statistically significant Fe emission. We 
also tested the sensitivity of the best-fit parameters to the 
energy ranges used for fitting the PN data by varying the 
low-energy boundary of the fits from 0.2 keV to 0.4 keV. 
We did not find any significant change in the parameters 
of the high-energy absorption lines. The exclusion of the 
0.2-0.4 keV PN data, however, did result in slightly poorer 
constraints on the ionization parameter of the gas caus- 
ing the low-energy absorption. A detailed description of 
the current calibration status of the XMM-Newton instru- 
ments can be found on the XMM-Newton SOC WWW 
site. A recent energy- dependent re-working of the PN re- 
sponse has resulted in a significant improvement in the 
calibration down to 0.2 keV and brings the PN and MOS 
detectors into better agreement. 

In Figure 2a we show the 68% and 90% confidence con- 
tours (based on model 3 of Table 2) for the photon indices 
versus column densities for the three epochs. The col- 
umn densities for the low-energy absorbers in models 3 
and 4 are consistent with no variation. A possible weak 
variation of the photon index between epochs 1 and 3 is 
detected. 68% and 90% confidence contours for the ion- 
ization parameter of the low-energy absorber versus col- 
umn density are presented in Figure 2b. The absorber 
in PG 1115+080 appears to have been more ionized dur- 
ing epoch 1. In particular, during epoch 1 the ionized 
absorber is characterized by an ionization parameter of 
£ = L/nr 2 = 145^ 2 er S cm s_1 and a hydrogen column 
density of A H = 1.33lg;35 x 10 22 cm ~ 2 , and during epoch 
2 by £ = 56^23 erg cm s _1 and a hydrogen column density 
of A H = l-41±o;27 x 10 22 cm" 2 , where L is the integrated 
5 eV-300 keV incident luminosity, n is the electron num- 
ber density of the absorber, and r is the distance between 
the absorber and ionizing source. 

Figure 2 suggests that the ionization parameter de- 
creased in the order of epochs 1, 3, and 2. This trend 
in ionization parameter is evident in the PN data shown 
in Figure 1. For epoch 1, as shown in Figure la, the spec- 
trum at low energies ~ 0.2-0.4 keV appears less absorbed 
relative to epochs 2 and 3, as expected for a more highly 
ionized absorber where the opacity from low-z metals is 
smaller. Additional confirmation for the trend in ioniza- 
tion parameter is provided in §3.1. 

Our spectral analysis of PG 1115+080 implies the exis- 
tence of two distinct absorbers with very different ioniza- 
tion parameters. In particular, if we interpret the appar- 
ent detection of absorption lines at rest-frame energies of 
7.27±g;?o keV and 9.79±?;j^ keV as being due to Fe, the 
most conservative assignment (giving the lowest outflow- 
ing velocity) is to highly ionized Fe XXV which requires 
log£ ^3.5. In contrast, the significant absorption features 
below rest- frame energies of ~ 1.6 keV are best fit with an 
ionization parameter of the order of log £ = 2 (see Table 
2). A justification of the interpretation of the absorption 
lines is given in §3 and presented in more detail in Char- 
tas et al. (2002, 2003). To illustrate better the absorption 
lines that can arise at these two distinct ionization levels 
we show in Figure 3 the expected absorption spectrum of 
PG 1115+080 in the 1 — 10 A rest-frame range as derived 



with the photoionization code XSTAR assuming a slab- like 
stationary absorber with Ah = 9 x 10 22 cm -2 , log£ = 2, 
and log£ = 3.5. The assumed column density of 9 x 
10 22 cm -2 lies between the estimated values of the column 
densities of the high ionzation absorbers of PG 1115+080 
listed in Table 4. At ionization levels of log£ ~ 3.5 and 
rest- frame wavelengths in the range 1.5-2 A (see panel a 
of Figure 3) the Fe lines are quite isolated and are there- 
fore apparent even with low-resolution spectroscopy. At 
rcst-framc wavelengths in the range of 6.2-10 A and at 
ionization parameters of log£ ~ 2 the spectra contain a 
larger density of strong absorption lines that arc difficult 
to resolve and identify (see panel b of Figure 3). 

The derived column densities for the low-ionization ab- 
sorber of PG 1115+080 (see Table 2) for the three epochs 
described in this paper are lower than those found in X- 
ray observations of typical BAL quasars. For example, 
the column densities of BAL quasars as inferred from X- 
ray observations lie in the range of 0.4-250 x 10 23 cm~ 2 
(e.g., Gallagher ct al. 2002 and Punsly 2006 present 
BAL-quasar column densities derived from X-ray obser- 
vations). We emphasize that most current column-density 
estimates derived from X-ray observations are poorly con- 
strained due to the limited S/N of most existing X-ray 
spectra of BAL and mini-BAL quasars. The presence of 
ionized or partially covering absorption cannot be accu- 
rately inferred from the analysis of individual X-ray spec- 
tra of most BAL quasars and if present would in gen- 
eral lead to an underestimate of the column densities. 
On the other hand, the large lensing magnifications of 
APM 08279+5255 and PG 1115+080 have provided mod- 
erate S/N X-ray spectra for these objects resulting in im- 
proved constraints compared to unlcnscd quasars. We note 
that the column density of PG 1115+080 appears to have 
varied significantly since an earlier observation. Specifi- 
cally, during a ROSATPSPC observation of PG 1115+080 
on 21 November 1991 the 0.2-2 keV X-ray flux of this ob- 
ject was found to be at a minimum level (low-state) of 
about 0.4 x 10~ 13 erg s _1 cm~ 2 (see Figure 6 of Char- 
tas 2000) and the intrinsic column was constrained to be 
1.2 ± 1.1 x 10 23 cm~ 2 . This column density is consis- 
tent with the typical values observed in BAL quasars. We 
estimated the optical-to-X-ray spectral slope, quantified 
as a ox = log(/ 2 kev// 2500 ^)/log(^ 2 kev/^ 250Q ^) (Tanan- 
baum, et al. 1979), where f 2 keV and / ° are the flux 

zoUUA 

densities at 2 keV and 2500 A in the quasar rest-frame, 
respectively. For the cases of unabsorbed(absorbcd) flux 
densities at 2 keV for the epoch 1 XMM-Newton observa- 
tion of PG 1115+080 we find a ox = -1.63(-1.67). We note 
that from the absorbed flux density at 2 keV for the 1991 
ROSAT PSPC observation (low-state) of PG 1115+080 
we find a ox — -1.97. Due to the limited quality of the 
ROSAT PSPC spectrum useful constraints on the unab- 
sorbed value of a ox could not be obtained for the 1991 
observation. For comparison the value of a ox expected for 
a typical AGN with a 2500 A luminosity density similar 
to that of the epoch 1 observation of PG 1115+080 based 
on the recent empirical relation of Steffen et al. 2006 (see 
equation 2 of Steffen et al.) is a ox — -1.58. The rms of a ox 
at log(L °) = 30.8 is about 0.2 (see Table 5 of Steffen 

2500 A 

et al.). We conclude that PG 1115+080 is slightly X-ray 
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weak during the three observation epochs presented in this 
paper but still within the scatter of a ox detected in typical 
AGN. 

3. DISCUSSION 

Our analyses of the Chandra and XMM-Newton spectra 
of PG 1115+080 and APM 08279+5255 presented here and 
in the discovery papers of Chartas et al. (2002, 2003) in- 
dicate that the high-energy broad absorption features de- 
tected in these objects are significant at the > 99.9% confi- 
dence level. Our claimed detections of relativistic winds in 
the X-ray spectra of PG 1115+080 and APM 08279+5255 
were based on the following steps: 

a) We chose a standard approach to fitting the X-ray 
spectra of PG 1115+080 and APM 08279+5255. In par- 
ticular, we first started with a simple model and added 
components of increasing complexity motivated both by 
the residuals of the spectral fits and our understanding of 
the physics and structure of mini-BAL and BAL quasars. 
The first model consisted of a simple power-law modified 
by intrinsic absorption. This choice was inspired by nu- 
merous empirical and theoretical studies that indicate that 
the main X-ray continuum emission component of quasars 
(e.g., rest- frame 2-30 keV) is produced by inverse Comp- 
ton scattering of soft photons from the accretion disk by 
hot electrons in a corona (e.g., Haardt & Maraschi 1991; 
Reeves & Turner 2000). An inspection of the residuals 
to fits with simple absorbed power-law models indicated 
the presence of significant complex absorption residuals at 
rest-frame energies above ~6.4 keV and rest-frame ener- 
gies below ~1.5 keV. 

b) To account for these fit residuals we considered a 
variety of spectral models. Several of these models were 
rejected in the discovery papers and this work. The list 
of rejected models includes the following: Broken power- 
law models with and without a soft blackbody component, 
models that included one and two absorption edges to ac- 
count for the high-energy absorption, models that included 
neutral absorption at the source to account for the low- 
energy absorption, and absorption by possible intervening 
material to account for the high-energy absorption. Single 
and double partial-covering models that try to fit the high- 
energy features of PG 1115+080 with iron edges require a 
curvature of the low-energy continuum that is inconsistent 
with our data. Additional data would be required to rule 
out even more complex partial-covering models. 

c) Motivated by the structure of the low and high- 
energy residuals, the fact that PG 1115+080 and 
APM 08279+5255 exhibit broad absorption lines in the 
UV (e.g., the properties of the UV broad absorption lines 
in PG 1115+080 and APM 08279+5255 have been pre- 
sented in Michalitsianos ct al. (1996) and Irwin et al. 
(1998), respectively), and previous claims of detections 
of low-energy ionized absorption in the X-ray spectra of 
BAL quasars, we fitted the spectra with a model that 
consisted of a power-law modified by ionized absorption 
at the source and two iron lines at the source to ac- 
count for the high-energy absorption. This model pro- 
duced acceptable fits to the X-ray spectra of PG 1115+080 
and APM 08279+5255 and resulted in a significant im- 
provement over all other models considered. We note 
that models that included high-energy absorption lines in 



PG 1115+080 and APM 08279+5255 resulted in improve- 
ments in the fits compared to models that did not include 
absorption lines at the greater than 99.9% confidence level. 
In the present work we also confirm this significant detec- 
tion by re-analyzing the same data with updated instru- 
mental calibration files. 

d) We made a plausibility argument based on the en- 
ergies of the identified lines, and the known energies of 
absorption lines from all the abundant elements, that the 
detected absorption lines are associated with highly ion- 
ized Fe Ka absorption. Based on this identification we 
estimated the velocities of the outflowing X-ray absorbing 
material. 

Several of the key derivations that follow in this section 
assume our interpretation that the high-energy absorption 
is due to lines arising from highly ionized Fe XXV and/or 
Fe XXVI and that the outflow velocities of these absorbers 
range between 0.1c and 0.4c. Here we provide a sum- 
mary of our justification of this interpretation that was 
presented in more detail in Chartas et al. (2002, 2003). 

Our analyses of the Chandra spectrum of 
APM 08279+5255 and the XMM-Newton spectrum of 
PG 1115+080 showed strong evidence for the presence of 
absorption lines at rest-frame energies of 8.05+ ' 08 keV and 
9.79±o ig keV for APM 08279+5255 and rest-frame ener- 
gies of 7.27ig'?S keV and 9.79^ ^ keV for PG 1115+080. 
For APM 08279+5255 we assumed that the two lines 
are produced by resonant absorption due to Fe xxv 
ls-2p (Fe XXVI Lya), and inferred that the X-ray ab- 
sorbers are outflowing with velocities of ~ 0.20c(0.15c) 
and - 0.4c(0.36c), respectively. For PG 1115+080 we 
assumed that the high energy absorption is due to two 
lines produced by resonant absorption due to Fe xxv 
ls-2p (Fe XXVI Lya), and inferred that the X-ray ab- 
sorbers are outflowing with velocities of ~ 0.1c(0.05c) and 
~ 0.4c(0.36c), respectively. The wind geometry that we 
assume is based on the unified BAL model (e.g., Wey- 
mann et al. 1991; Murray et al. 1995; Proga 2000; Elvis 
2000). In this model the AGN outflow originates in the 
accretion disk from a narrow range of radii. This wind 
rises initially almost perpendicular to the accretion disk 
and becomes more radial and equatorial at larger radii 
to form a bi-cone. The unified BAL model proposes that 
most of the observed range of absorption line-widths can 
be explained with orientation and with a velocity gradient 
in the outflowing stream. The hydrodynamical models of 
Proga et al. (2000) indicate that the density of the wind 
peaks near an inclination angle of about 70 degrees and 
most of the outflow is confined within ± 10 degrees of this 
angle. We note that Proga & Kallman (2004) showed that 
the wind inclination angle can be smaller and the wind 
does not need to be equatorial if it is mostly driven by UV 
radiation from the accretion disk. 

For the velocity calculations we considered the special- 
relativistic velocity correction and assumed that the angle 
9 between the wind velocity and our line of sight is 20°. 
This angle is not constrained with the present data; how- 
ever, hydrodynamical simulations indicate that the BAL 
wind divergence angle may range between 10°-30° depend- 
ing on the location of the inner radius of the disk. The 
velocity estimates are not sensitive to this range of an- 
gles. For example, the inferred outflow velocity (assum- 
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ing 9 ~ 20°) for the 9.79 keV absorption line detected in 
PG 1115+080 would vary from -6% to +13% for 6 ranging 
between 10°-30°. 

Of all the abundant elements, iron absorption lines 
would be the closest in energy to the observed features. 
As shown in Figure 3 the strongest lines with rest-frame 
energies indicated in parenthesis near the observed ones 
are: Fe XXV ls-2p (6.7 keV), Fe XXV ls-3p (7.88 keV), 
Fe XXVI Lya (6.97 keV), Fe XXVI ls-3p (8.25 keV), S 
XV ls-2p (2.46 keV), S XVI Lya (2.62 keV) Si XIV Lya 
(2.005 keV), Ar XVIII Lya (3.321 keV), and Ca XX Lya 
(4.104 keV) (based on the energies of permitted resonance 
lines of Verner et al. 1996). In this sense, our interpreta- 
tion that the absorption lines are associated with highly 
ionized Fe K absorption is the most conservative one pos- 
sible (e.g., absorption lines from relativistic sulfur or oxy- 
gen would require much larger blucshifts). We also inves- 
tigated whether intervening absorbers in the lens galax- 
ies or in possible damped Lyman alpha systems along the 
line of sight could explain the high-energy absorption fea- 
tures. We concluded that the observed high-energy ab- 
sorption features cannot be produced by absorption in 
intervening systems and the most likely origin is intrin- 
sic absorption by highly ionized iron. The intrinsic origin 
of the absorption has been confirmed with the detection 
of significant flux variability of the absorption features in 
APM 08279+5255 (see Figure 7 in Chartas et al. 2003) 
and now with the detection of significant variability in 
PG 1115+080 (see Figure 4) as discussed in more detail in 
§3.1. From the present data we can infer significant vari- 
ability of the normalization of the absorption lines, how- 
ever, higher quality spectra will be required to distinguish 
if there is any energy variability of the absorption lines. 

Assuming this interpretation for the relativistic X-ray 
absorbing material, we present in the following sections 
estimates of the mass-outflow rate and efficiency of the 
outflow in PG 1115+080 and APM 08279+5255. These 
estimates are important for addressing the fundamental 
issue of whether quasar outflows are capable of influenc- 
ing their host galaxies and black-hole growth. Since the 
present observations of PG 1115+080 do not resolve the 
lensed images we also estimate the influence of this effect 
on our main conclusions. 

3.1. Variability of the outflow in PG 1115+080 

Our spectral analysis of the XMM-Newton observations 
of PG 1115+080 presented in §2.1 indicates possible vari- 
ability of the X-ray broad absorption features between 
different epochs. To quantify the significance of these 
variations we determined the sum of the residuals of Ax 
between the best-fit 2-2.5 keV and 5-10 keV continuum 
model and data in the energy range 2.5-5.0 keV based on 
the spectral fits presented in Figure 1. As described in §2.2 
we fit the spectra within the energy ranges of 2-2.5 keV 
and 5-10 keV with Galactic absorption and a power-law 
model and extrapolated this model to the energy ranges 
not fit (see §2.1). For this summed-Ax method to be 
statistically valid, the best-fit 2-2.5 keV and 5-10 keV 
models and the S/N of the spectra need to be similar for 



all epochs. The best-fit 2-2.5 keV and 5-10 keV spectral 
slopes for all three epochs have essentially the same value 
of r = 1.9 for epochs 1, 2, and 3 respectively, and the S/N 
of the three observations are similar with effective expo- 
sure times and source count rates listed in Table 1. We 
therefore do not expect variations of the best-fit models or 
S /N of the spectra to contribute strongly to the variations 
in the high-energy residuals. We find total A% residuals 
of A X Ei = -40 ± 3, A XE2 = -26 ± 4, and A XE3 = -20 
± 3 for epochs 1,2, and 3, respectively. 

A second approach to quantifying variability of the high- 
energy absorption features relies on taking ratios of the 
spectra. This approach is independent of model assump- 
tions and takes into account the uncertainties in each bin. 
For the purpose of this analysis the grouping of the spectra 
was made identical for the three epochs. In Figure 4 we 
show the ratios of the spectra of epoch 1 to epoch 2 (R12), 
epoch 1 to epoch 3 (-R13), and epoch 2 to epoch 3 (i?23)- 
We find the mean of the ratios within the high-energy ab- 
sorbed range 2.5-5.0 keV to be < R 12 > = 0.83 ± 0.05, 
< R 13 > = 0.73 ± 0.05, and < i? 23 > = 0.88 ± 0.05. 

As a third method of estimating the variability of the X- 
ray broad absorption features we compared the strengths 
of the best-fit Gaussian lines for the three epochs. Specif- 
ically, we found the total absorbed photon fluxes in the 
Gaussian absorption lines for epochs 1, 2, and 3 were 
(2.4 ± 1.1) x 10~ 5 photons cm~ 2 s -1 , (0.49 ± 0.20) x 
10~ 6 photons cm~ 2 s -1 , and (0.34 ± 0.16) x 10~ 6 pho- 
tons cm -2 s _1 , respectively. The decrease in the total ab- 
sorbed photon fluxes in the absorption lines is consistent 
with the results of the other two methods. 

All approaches imply significant variability of the X-ray 
broad absorption features in PG 1115+080 between epochs 
1 and 2 separated by 0.92 yr (rest-frame) and marginal 
variability of the X-ray broad absorption features between 
epochs 2 and 3 separated by 5.9 d (rest-frame). The 5.9 d 
variability, if real, is consistent with a relatively small 
launching radius of the X-ray absorber. Specifically, as- 
suming a radiatively driven wind and using equation 1 
of Chartas et al. (2003) we estimate that a launching ra- 
dius of about 7 R s is needed for the absorber to reach a 
terminal velocity of about 0.4 c. At a launching radius 
of 7 Rg the time to reach 90% of the terminal velocity is 
about 9 d. This time-scale is consistent with the marginal 
detection of variability of the Fe absorption lines between 
epochs 2 and 3. 

A model-independent confirmation of the trend of the 
ionization parameter presented in §2.1 can also be seen in 
Figure 4. In particular, Figure 4 shows significant positive 
residuals of the ratios R12 and R13 and negative residuals 
of the ratio R23 within the 0.2-0.4 keV band. This indi- 
cates that the 0.2-0.4 keV flux decreased in the order of 
epochs 1,3, and 2. This trend in soft X-ray flux is consis- 
tent with the observed trends of the 0.2-0.4 keV residuals 
shown in Figure 1 and the observed trend in ionization 
parameter derived from our spectral analysis and shown 
in Figure 2. We note that the ratio of the EPIC PN spec- 
tra is not sensitive to the precise value of the calibrated 
effective area (the effective area cancels out in the ratio) 

5 We find the bolometric luminosity of PG 1115+080 to be Lb i = 3.3 X 10 46 erg s _1 . To estimate Lb i we used the empirical relation Lb i 
= /bc(3000A)^ L„(3000A)/j -1 , where /bc(3000A) = 5.3 is the luminosity-dependent bolometric correction obtained from equation (21) of 
Marconi et al. (2004) and fi is the flux magnification assumed to be ~ 25 based on lensing models of this system (e.g., Impey et al. 1998). 
The value of L„(3000A) = 1.57 X 10 32 erg s" 1 Hz" 1 for PG 1115+080 was taken from Neugebauer et al. (1987). 
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as long as the response of the PN detector remained the 
same between observations. The PN detector response is 
known to be stable with time according to the calibration 
status of the XMM-Newton instruments (see footnote 4) . 

Significant variability of PG 1115+080 has also been re- 
ported in the optical, UV and X-ray bands. Wcymann 
et al. (1980) report several absorption lines as being vis- 
ible in the May 1980 observations of this quasar. Young 
et al. (1982) detect only weak absorption features blue- 
wards of the CIV line in their December 1980 spectrum of 
PG 1115+080 and do not detect the lines seen in the May 
1980 observations. Strong and rapid variability of broad 
absorption lines blueward of OVI in PG 1115+080 were 
reported by Michalitsianos et al. (1996). Chartas et al. 
(2000) reported a decrease by a factor of about 13 of the 
0.2-2 keV flux between 1979 December 5 and 1991 Novem- 
ber 21 and an increase by a factor of about 5 between the 
1991 November 21 and 1994 May 27 observations of PG 
1115+080. 

Several recent theoretical studies of high-velocity out- 
flows from quasars (e.g., Everett & Ballantyne 2004; Sim 
et al. 2005) assume that the high-ionization absorption 
lines of S XIV, S XVI and Fe XXV observed in several 
quasars (i.e., PG 1211+143 and PG 0844+349 reported 
by Pounds at al. 2003a, b) originate from gas that has 
reached its terminal velocity. This terminal velocity is 
then approximated with the escape velocity from the re- 
gion from which the wind is launched, resulting in the 
approximation i?i aU nch ~ Rs{c/v b a ) , where w b s is the 
observed outflow velocity and R s is the Schwarzschild ra- 
dius. For u bs in the range 0.1-0. 4c the expected launching 
radii for PG 1115+080 will lie in the range 100 - 6i? s . For 
an estimate of the black-hole mass we used the empirical 
relation (equation 7) of Vestergaard & Peterson (2006). 
Specifically, we find the mass of the black hole and the 
Schwarzschild radius to be Mbh = 1.1 x 10 9 M Q and R s 
= 2GMbh/c 2 = 3.1 x 10 14 cm, respectively, based on the 
FWHM of CIV of 4700 km s" 1 taken from the Hale spec- 
trum of PG 1115+080 (Young et al. 1982), and the 1350 A 
luminosity density of L v = 9.48 x 10 31 fi^ 1 erg s _1 Hz -1 
taken from Neugebauer et al. (1987), where fi is the flux 
magnification assumed to be ~ 25 based on lensing models 
of this system (e.g., Impey et al. 1998). 

One of the assumptions in these recent theoretical stud- 
ies is that the absorbers are observed near their termi- 
nal velocities. As we proposed in our previous study of 
APM 08279+5255 another possibility is that we are ob- 
serving the absorber as it is being accelerated near the 
launching radius. We speculate that initially the absorb- 
ing material supplied by the accretion disk has a relatively 
low ionization parameter because of the large gas density 
at the base of the wind. A low ionization parameter re- 
sults in a large value of the force multiplier that can lead to 
significant acceleration through scattering in atomic reso- 
nance lines (line-driving; the force multiplier represents the 
ratio by which the line opacity of the absorbing material 
increases the radiation force relative to that produced by 
Thomson scattering alone). Because of the intense radia- 
tion field of the central source the absorbing gas becomes 
increasingly ionized within a time that is much shorter 
than the time it takes for the absorber to reach its termi- 
nal velocity. When the ionization parameter reaches values 



of about £ > 1000 erg cm s _1 absorption from highly ion- 
ized lines such as Si XIV, S XVI and Fe XXV are expected 
to be produced. As the ionization parameter increases the 
absorber becomes completely ionized and line-driving be- 
comes no longer important. This interpretation predicts 
rapid changes (compared to those observed in UV broad 
absorption lines) in the observed energies and equivalent 
widths of the X-ray BALs. We note that magnetically 
driven outflows is also a possibility, as this might be needed 
to reach such high velocities if the material becomes highly 
ionized quickly. 

The average continuum flux of PG 1115+080 has not 
varied by more than ~ 20% between epochs as indicated 
by the observed PN count rates listed in Table 1. The 
fractional change of the continuum flux can also be seen 
in the ratio plots shown in Figure 4. There is thus no 
correlation between the depth of the high-energy broad 
absorption features and continuum flux level. 

3.2. Efficiency of the quasar outflow in PG 1115+080. 

One of the key and highly uncertain parameters used in 
recent theoretical models that describe black-hole growth 
and structure formation is the efficiency of the quasar out- 
flow. The efficiency is defined as the fraction of the total 
bolometric energy released over a quasar's lifetime into the 
ISM and IGM in the form of kinetic energy injection and 
can be expressed as: 

e k ,i = = 2nf c , i R i (R i /AR i )N^m p ^ (1) 

^ -^Bol -^Bol 

where M is the mass-outflow rate of component i, w w ind,i 
is the outflow velocity of the X-ray absorber of compo- 
nent i, / Cj j is the global covering fraction of the absorber 
of component i, ARi is the thickness of the absorber at 
radius R[ of component i, Nu,i is the hydrogen column 
density of component i, and LboI is the bolometric pho- 
ton luminosity of the quasar. We focus on estimating ek 
for the outflow in PG 1115+080 observed during epoch 1. 
The velocities of the different outflowing X-ray absorbers 
are listed in Table 3. We assume a conservatively wide 
range for the covering factor of / c =0. 1-0.3 based on the 
observed fraction of BAL quasars (e.g., Hewett & Foltz 
2003) and a fraction R/AR ranging from 1-10 based on 
current theoretical models of quasar outflows (e.g., Proga 
& Kallman 2000). As we argued in § 3.1 we are likely 
observing the absorbers as they are accelerated near their 
launching radii. We have assumed absorber radii ranging 
from 3-15 R s based on the maximum observed X-ray ve- 
locity components. We used a Monte Carlo approach to 
estimate the errors of M and £k- The values of f c , R/AR, 
and Ri were assumed to have uniform distributions within 
their error limits. By multiplying these three distribu- 
tions and with the appropriate constants from equation 1 
we obtained the distributions of Mi and e^. We finally de- 
termined the means of the distributions of M and Ck and 
estimated the 68% confidence ranges. 

We note that the last stable orbit can be even smaller 
in a Kerr black hole and the disk can extend down to 
the event horizon. This may also be possible even for 
a Schwarzschild black hole, as argued recently by Agol 
& Krolik (2000). However, at radii smaller than 3 R s 
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the gravitational redshift would begin to become impor- 
tant with E b s /Ee m is ~ - 2M b h/r C mis) and one would 
have to invoke even larger outflow velocities to explain the 
detected blue-shifted Fe lines. 

The total hydrogen column densities of the Fe XXV Ka 
components absl and abs2 are consistent with the ones de- 
termined in Chartas et al. (2003). The mass-outflow rates 
and efficiencies of the outflowing components are listed in 
Table 4. 

There are several systematic errors in estimating the 
radii and column densities of the absorbers that are diffi- 
cult to constrain and have not been included in the present 
calculations of the outflow efficiencies. However, both 
these systematic uncertainties are one sided in the sense 
that they lead to underestimates of the outflow efficiencies. 
Specifically, we have allowed the radii of the absorbers to 
range between the last stable orbit of the black hole to 15 
R s based on the maximum inferred outflow velocities. We 
therefore expect any errors in the estimates of the radii of 
the absorbers to lead to larger values of outflow efficien- 
cies. Possible saturation of the absorption lines will lead 
to an underestimate of the column densities and therefore 
to an underestimate of the outflow efficiencies as well. 

3.3. Contamination from unresolved lensed images 

The lensed images of PG 1115+080 cannot be resolved 
with XMM-Newton, and therefore the spectra used in our 
analyses contain contributions from all images. We esti- 
mate how this effect influences our previous results by first 
determining X-ray flux ratios based on recent Chandra ob- 
servations of PG 1115+080 that do resolve the lensed im- 
ages. To estimate the X-ray flux ratios of PG 1115+080 
we modeled the Chandra images of Al, A2, B and C, with 
point-spread functions (PSFs) generated by the simulation 
tool MARX (Wise et al 1997). The X-ray event locations 
were binned with a bin-size of 0"0246 to sample the PSF 
sufficiently (an ACIS pixel subtends 0"491). The simu- 
lated PSFs were fitted to the Chandra data by minimizing 
the C-statistic formed between the observed and simulated 
images of PG 1115+080. The relative positions of the 
images were fixed to the observed NICMOS values taken 
from the CfA- Arizona Space Telescope Lens Survey (CAS- 
TLES). In Figure 5 we show the Lucy- Richardson decon- 
volved image of the Chandra observation of PG 1115+080. 
We find that the relative X-ray flux ratios in the 0.2-8 keV 
band are [Al/C] fuU = 3.86 ± 0.42, [A2/C] full = 1.85+ 0.22, 
and [B/C]fuii = 0.93 ± 0.14. For comparison, the HST H- 
band flux ratios are [A1/C] H = 4.06 ± 0.17, [A2/C] H = 
2.56 ± 0.12, and, [B/C] H = 0.65 ± 0.04. The #-band 
magnitudes were taken from the CASTLES WWW sited 
The X-ray flux ratios of PG 1115+080 indicate that about 
75% of the X-ray flux originates from the close-separation 
images Al and A2. Since the observed time delay between 
images Al and A2 is Ai A iA2 = 0.149 ± 0.006 d there may 
be variability of the X-ray BALs within this time-scale 
only if the absorbers are launched near the last stable or- 
bit of the black hole. The "long" time delays between 
images were presented in Schechter et al. (1997) and have 
image C leading images A1+A2 by 9.4 ± 3.4 d and image 
C leading image B by 23.7 ± 3.4 d. Since images B and 
C are considerably weaker than A1+A2 they are unlikely 

6 The CASTLES WWW site is located at |http: //cfa-www.harvarcie 



to be able to produce a significant change in the combined 
spectrum of PG 1115+080. 

3.4. Comparison with APM 08279+5255 

BAL quasars are relatively weak in X-rays making the 
detection of broad absorption features difficult with XMM- 
Newton and Chandra. We have previously presented evi- 
dence for a relativistic outflow in the gravitationally lensed 
z=3.91 BAL quasar APM 08279+5255. Following an ap- 
proach similar to the one described in §3.1 and using 
the derived outflow properties presented in Chartas et 
al. (2003), we find that the mass-outflow rates and ef- 
ficiencies of the outflows associated with the two Fe ab- 
sorbers in APM 08279+5255 are M abs i = 1.7±^ M 

yr" 1 , e k ,absi = 10±| x 10- 3 , M abs2 = 3.3±£? M yr- 1 , 
and £k,abs2 = 8^5 x 10~ 2 . The total efficiency of the 
outflow will be the sum of all components of all species 
over time. Because of the low S/N of the available spec- 
tra of APM 08279+5255 additional components cannot be 
constrained, and we should consider the estimated total 
mass-outflow rates and efficiencies as lower limits. The es- 
timated outflow rate is an order of magnitude smaller than 
the accretion rate in APM 08279+5255 which we estimate 
to be 1.8 x l0- 3 (L 44 /n)M Q yr" 1 - 40 M Q yr" 1 , where 
we assumed a typical accretion efficiency of r] = 0.1. 

4. CONCLUSIONS 

High-energy X-ray absorption is detected in three 
epochs of monitoring observations of the mini-BAL quasar 
PG 1115+080. We interpret this absorption as due to 
lines arising from highly blueshifted Fe XXV ls-2p and/or 
Fe XXVI Lya. We find that the depths of the X-ray 
broad absorption features decreased significantly between 
epochs 1 and 2 separated by 0.92 yr (proper-time) and 
detect a marginal decrease between epochs 2 and 3 sepa- 
rated by 5.9 days (proper-time). The case for a relativis- 
tic flow on the basis of either the epoch 2 or 3 observa- 
tions of PG 1115+080 alone is less compelling. We note 
that models that included high-energy absorption lines 
in PG 1115+080 and APM 08279+5255 resulted in im- 
provements in the fits compared to models that did not 
include absorption lines at the greater than 99.9% confi- 
dence level. We had previously reported rapid variability 
over timescales of 1.8 weeks (proper-time) of X-ray BALs 
in the quasar APM 08279+5255. Variability of similar 
magnitude and over similar time-scales has never been ob- 
served in UV BALs. 

Assuming the interpretation that the high-energy ab- 
sorption features are due to lines arising from highly ion- 
ized Fe XXV we used the measured outflow velocities, 
column densities, and estimated launching radii to con- 
strain the mass-outflow rates and outflow efficiencies for 
PG 1115+080 and APM 08279+5255. We find the out- 
flow efficiencies to be £k = 0.64^0^(68% confidence), and 
ek = 0.09^05(68% confidence), respectively. These esti- 
mates include only contributions from observed compo- 
nents and therefore should be considered as lower lim- 
its. Our derived estimates of the efficiency of the out- 
flows in mini-BAL quasar PG 1115+080 and BAL quasar 
APM 08279+5255, when compared to values predicted by 

ensdata/ 
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recent models of structure formation (SO04; G04; SDH05), 
imply that these winds will have a significant impact on 
shaping the evolution of their host galaxies and in regu- 
lating the growth of the central black hole. 

We note that our reported values of the mass-outflow 
rates and outflow efficiencies represent instantaneous 
quantities. Because of the detected variability of the out- 
flows in PG 1115+080 and APM 08279+5255 it is appro- 
priate to compare the outflow efficiency incorporated in 
theoretical models with the average of this observed quan- 
tity over a period that covers the range of variability. Ad- 
ditional monitoring observations are needed to constrain 
better the variability of the outflow properties. 

The fraction of BAL quasars with highly blucshiftcd 
Fe absorption lines is not well constrained from avail- 
able X-ray observations because of the relatively poor- 
to-moderate S/N X-ray spectra of BAL quasars ob- 
served to date. We have shown that the two X-ray 
brightest mini-BAL and BAL quasars, APM 08279+5255 
and PG 1115+080, show significant highly ionized and 
blucshifted Fe lines. The third brightest BAL quasar 
known to date CSO 755 shows a hint of an Fe absorption 
feature in the rest-frame energy range 7.1-9.3 keV (Shem- 
mer et al. 2005). Higher quality observations are required 
to confirm this result. The X-ray spectra of the remaining 
observed BAL quasars are of lower S/N and do not pro- 
vide interesting constraints on the presence of high-energy 
Fe absorption features. 

To date, the most interesting constraints on quasar out- 
flows are derived from X-ray and UV spectroscopic anal- 
yses. Current observations of BAL quasars imply that 
the X-ray and UV bands may be sampling different parts 
of the absorbing outflows. A large fraction of the high 
velocity BAL outflow is apparently visible only in the X- 
ray band. Unfortunately, the limited number and qual- 
ity of available X-ray spectra leave large uncertainties in 
the mass-outflow rate. Absorption studies of quasars such 
as PG 1115+080 and APM 08279+5255 with future high 
resolution, high throughput X-ray spectroscopic missions 
such as Constellation- X hold great promise for pinning 
down this crucial quantity. 
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Table 1 

Log of Observations of BAL Quasar PG 1115+080 









Exposure 


Effective 




Observation 


Observatory 


Observation 


Time 


Exposure Timc a 




Date 




ID 


(ks) 


(ks) 




2004 June 26 


XMM-Newton 


0145750101 


86.52 


54.92 


0.181 ± 0.002 


2004 June 10 


XMM-Newton 


0203560201 


86.65 


51.46 


0.212 ± 0.002 


2001 November 25 


XMM-Newton 


0082340101 


63.36 


52.37 


0.202 ± 0.002 


2000 November 3 


Chandra 


1630 


9.95 


9.83 


0.084 ± 0.003 


2000 June 2 


Chandra 


362 


26.83 


26.5 


0.075 ± 0.002 



"Effective exposure time is the time remaining after the application of good time-interval (GTI) tables to remove portions of the observation 
that were severely contaminated by background. 

'Background-subtracted source count rate including events with energies within the 0.2-10 keV band. The source count rates and effective 
exposure times for the XMM-Newton observations refer to those obtained with the EPIC PN instrument. Please see §2 for details on source 
and background extraction regions used for measuring R sc . 
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Table 2 

RESULTS FROM FITS TO THE XMM-Newton EPIC SPECTRA OF PG 1115+080 



Model" Parameter Values For Epoch l c Values For Epoch 2 C Values For Epoch 3 C 

2001 November 25 2004 June 10 2004 June 26 

1 r 1.85ir 3 3 l-85±°o°oI l-80iS:°o 3 3 

AT H (0.221°;°*) X 10 22 cm- 2 (0.461°;°b) X 10 22 cm" 2 (0.37l°;°f) X 10 22 cm" 

X 2 /v 361.2/240 361.1/272 344.8/255 

P( x 2 /v) d 6.6 x 10~ 7 2.3 x 10~ 4 1.5 x 10~ 4 

9 p i co + 0-04 , no + 0.03 , 79 + 0.03 

1 - b8 -0.04 i-^-o.oa 1 - 72 -0.03 



Nn (0.12±°-_°|) x 10 22 cm- 2 (0.44l°,;°,|) X 10 22 cm" 2 (0.32l°,°|) x 10 2 

E absl 7.271°, H keV 6A0±°\* keV 

cr abB i < 0.50 keV < 0.44 keV 

EW ab3l 0.1171°,;°?? keV 0.0571°,;°,™ keV 

E aba2 10.5C" ™ keV 8.53i° ; 3 ° keV 9-6±o. 78 kcV 

<x abB2 3.49± 4 ;|| keV < 0.53 2.7l 2 ; 7 keV 

EW aba2 2.271°+? kcV 0-!8-o 10 kcV 1 -^-ol k °V 

X 2 /v 284.8/228 320.8/260 363.6/250 

P(X 2 Ivf 6.27 x 10~ 3 6.04 x 10~ 3 3.47 x 10~ 6 



1 i - 84 -0.03 1 -"-0.04 i-88-0.04 

at 1 1 qq + 0.58^ w i n 22 _-2 /-, ,..|+0.36\ „ i n 22 -2 /-, r; o + 0.46* w i n 22 _-2 

Wh (1.33_ 35 ) X 10 cm (1.41_ 27 ) X 10 cm (1.52_ 37 ) x 10 cm 

£ 144.9 + 5 3 3 Cr S cm s_1 55 - 9 t23 4 Cr S Cm S_1 90 - 2 i38 5 Cr S Cm S_1 

E ab3l 7.27+° ™ keV 6.40+° ;! keV 

cr abB i < 0.44 keV < 0.33 kcV 

EW aba i 0.1131°;°?? keV 0.0691°;°™ kcV 

E aba2 9.79l° °° keV 8.6ll°; 2 ? keV 8-91q.I1 kcV ( 68 % lcvcl ) 

cr abB2 3.07lJ | 2 keV < 0.51 keV O-SIq.L kcV ( 68 % level) 

EW aba2 2.05l° .4? kcV 0.17l° J keV 0.281°?? keV (68% level) 

x 2 /v 241.5/227 265.1/259 278.8/249 

P(x 2 /v) d 0.20 0.38 0.10 

T 1 84+ 05 1 94+ 02 1 90+°-° 4 

AT /i fV7 + - 70 ^ w ir>22 m -2 /-, C1 +0.3CK „ i n 22 -2 /-, « Q + 0.39\ w i n 22 _-2 

ivh (l-97_ 49 ) x 10 cm (1.61_ 30 ) x 10 cm (1.731 37 ) X 10 cm 

/ 0.40l°;° ° 0.65l° ;° 3 0.57l°;° 4 4 

E abal 7.27l° 2 | keV 6.41±g"i| kcV 

cr abB i < 0.44 keV < 0.33 keV 

EW aba i 0.103l° °? 2 keV 0.075l° °™ keV - 

E aba2 9.62+1^ kcV 8.53l°; 2 ° kcV s - 9 -p.| kcV ( 68 % level) 

3.181°?? kcV < 0.44 keV 0.5+ ~ 

llSII k <=V 0.191°, ■{{ kcV 0.28l 013 

X 2 /v 247.6/227 263.9/259 279.6/249 

P(x 2 jv) d 0.17 0.40 0.09 



"■ a bs2 

EW aba2 2.011°-^ keV 



a Model 1 consists of a power law and neutral absorption at the source. Model 2 consists of a power law, neutral absorption at the source, and Gaussian 
absorption lines at the source. Model 3 consists of a power law, ionized absorption at the source, and Gaussian absorption lines at the source. Model 
4 consists of a power law, partial covering of the source, and Gaussian absorption lines at the source. All model fits include the Galactic absorption 
toward the source (Stark et al. 1992). All fits have been performed on the combined spectrum of images Al, A2, B and C of PG 1115+080. 

All absorption-line parameters are calculated for the rest frame. 
C A11 errors are for 90% confidence unless mentioned otherwise with all parameters taken to be of interest except absolute normalization. 
d P(x 2 /v) is the probability of exceeding \ 2 f° r v degrees of freedom if the model is correct. 
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Table 3 

Energies and Inferred Outflow Velocities of Absorption Lines Detected in the XMM-Newton EPIC 

Observations of PG 1115+080 During Epoch 1 



LillC " Component £ ob s Brest Blab EWiost Croat »ab! 

(keV) (keV) (keV) (keV) (keV) (c) 

Fe XXV ls-2p absl 2.67±g;J* 7.27+°;^ 6.7 0.11^° ™ <0.44 0.09+°;^ 

Fe XXV ls-2p afes2 3.60lg;%g 9.79±° 1 9 6 5 6.7 2.05±g;|^ 3.07±^;|| 0.40±g;°° 



Note. — The present spectra cannot distinguish whether the absorption in components absl and abs2 is caused by Fe XXV ls-2p or 
Fe XXVI Lya. 
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Table 4 

Hydrogen Column Densities, Mass-Outflow Rates and Efficiencies of Outflows in PG 1115+080 and 

APM 08279+0552 



Component 


Lino N H M 

cm- 2 (M yr' 1 ) 

PG 1115+080 Outflow 






absl 
abs2 


Fo XXV ls-2p 4 x 10 22 0.101°°:! 
Fe XXV ls-2p 4 x 10 23 4.6+* ° 


7.2+« 3 x 
6.3lt 5 o >< 


io~ 4 
io- 1 


APM 08279+5255 Outflow 


absl 


Fo XXV ls-2p 1 x 10 23 l-7t\.i 


i.o±°;l x 


10~ 2 


abs2 


Fe XXV ls-2p 1 x 10 23 3.3± 2 ° 


0.8+0 5 x 


io- 1 



Note. — The bolometric luminosites of PG 1115+080 and APM 08279+0552 are 3.3 X 10 46 erg s" 1 and 2 x 10 47 erg s _1 , respectively. 
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Rest-Frame Energy (keV) Rest-Frame Energy (keV) 

0.544 1.36 2.72 5.44 13.6 27.2 0.544 1.36 2.72 5.44 13.6 27.2 




Observed-Frame Energy (keV) Observed-Frame Energy (keV) 



Fig. 1. — (a) The top panel shows the XMM-Newton PN spectrum of the combined images of PG 1115+080 for Epoch 1 fit with Galactic 
absorption and a power-law model to events with energies lying within the observed-frame ranges of 2-2.5 keV and 5-10 keV. The best-fit 
values of the power-law photon indices in these energy ranges for the three epochs were almost identical with T ~ 1.9. The lower panel shows 
the residuals of the fit in units of lc deviations. Several absorption features within the observed-frame range of 1.5-5.2 keV are noticeable in 
the residuals plot, (b) Same as (a) for Epoch 2. (c) Same as (a) for Epoch 3. For clarity we only show the higher S/N ratio PN data; however, 
all fits were performed simultaneously using the PN and MOS1+2 data unless mentioned otherwise. The vertical dashed lines indicate the 
best-fit energies of the high-energy absorption lines for epoch 1. The top panels (d), (e) , and (f) show the same data shown in panels (a), (b), 
and (c) overplotted with the best-fit models taken from model 3 of Table 2. The lower panels of (d), (e), and (f) show the reduced residuals 
of these fits. 
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Fig. 2. — 68% and 90% confidence contours between (a) the photon index and the intrinsic-absorption column density, and (b) the 
ionization parameter and the intrinsic-absorption column density for the XMM-Newton observations of PG 1115+080 during epoch 1 (blue 
dotted contours), epoch 2 (green dashed contours) and epoch 3 (black solid contours) assuming Model 3 of Table 2. The ionization parameter 
and column density are for the low-energy absorber. 
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Rest-Frame Energy (keV) 

6.198 3.099 2.066 1.549 

i 1 1 1 1 1 1 




o 

Rest-Frame Wavelength (A) 



Fig. 3. — Simulations of absorbed spectra with the continuum level normalized to that of PG 1115+080. In panels a and b we show the 
absorbed spectra for ionization parameters of log£= 3.5 and log £=2.0, respectively. 
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Observed-Frame Energy (keV) 



Observed-Frame Energy (keV) 



Fig. 4. — The ratio of the PG 1115+080 spectra of (a) epoch 1 to epoch 2, (b) epoch 1 to epoch 3 and, (c) epoch 2 to epoch 3. The spectra 
are binned with the same grouping used for epoch 1. The spectrum of epoch 1 was binned such that the minimum counts per bin was 100 
counts. Panels (d), (e), and (f) on the right show the same ratios as the panels on the left with the only difference being that the spectrum 
of epoch 1 was binned such that the minimum counts per bin was 200 counts. We note that the PN background for these observations 
becomes increasingly significant compared to the source signal for observed-frame energies above 6.5 keV. This explains the larger error bars 
and dispersion of the data points for the ratios above 6.5 keV. 
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1 arcsec 





Al 



Fig. 5.— 



Lucy-Richardson deconvolved image (0.2 - 8 keV band) of the Chandra observation of PG 1115+080. 



